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Abstract 

In supersymmetric models whose gauge group includes an addi- 
tional U{1) factor at the TeV scale, broken by the VEV of an standard 
model singlet S, the parameter space can accommodate a very light 
neutralino not ruled out experimentally. This higgsino-like fermion, 
stable if i?-parity is conserved, can make a good cold dark matter 
candidate. We examine the thermal relic density of this particle and 
discuss the prospects for its direct detection if it forms part of our 
galactic halo. 
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1 Introduction 



Supersymmetric models with conserved i?-parity have in the hghtest super- 
symmetric particle (LSP) a natural candidate for dark matter. This very 
appealing feature has motivated much work (see |^ for review and refer- 
ences) and many increasingly sophisticated studies on the field. Most of these 
analyses have concentrated in the minimal supersymmetric standard model 
(MSSM) with or without theoretical constraints on its wide parameter space. 
It is important, however, to keep an open mind to the possibility that this 
simplest supersymmetric generalization of the standard model may not be 
the one realized in nature. One should be careful not to identify the predic- 
tions of the MSSM with those generic of low-energy supersymmetry. In this 
respect it is healthy to explore (well motivated) extensions of the MSSM in 
search of phenomenological (or cosmological) consequences that are different 
from those expected in the minimal model. 

Perhaps among the best motivated extended models are those that include 
an additional U (1) factor in the gauge group, broken radiatively at some scale 
below 0{1) TeV by the VEV of a standard model singlet S. This type of 
models can be generically expected to arise as the low-energy limit of some 
string models and have a number of interesting consequences both at the 
phenomenological (e.g. for Z' and Higgs physics) and theoretical level (e.g. 
they can accommodate a natural solution to the /^-problem). For a detailed 
study of this kind of scenarios we refer the reader to refs. 0, In this 
letter we would like to consider the lightest neutralino in this type of models 
as a possible good dark matter candidate. We do not attempt a complete 
exploration of the parameter space of these models (even wider than that 
of the MSSM) but rather focus on a particular region in which the LSP 
has properties completely different from those that could be expected in the 
MSSM. As is described in the next section, in the region we study, the LSP 
is a neutralino mainly composed of the fermionic superpartner of the singlet 
used to break the extra f/(l). The relic density of such particle, computed in 
section 3, turns out to be of the right order of magnitude for an interesting 
dark matter candidate, and the prospects for its laboratory detection are 
estimated in section 4. Finally, in section 5 we present some conclusions. 
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2 The Dark Matter Candidate 



The symmetry breaking sector in these models includes a chiral multiplet S, 
singlet under the standard model gauge group but with a U{iy charge Qs- 
This field couples in the superpotential 



W = hSHiH2 



(1) 



to the usual Higgs doublets Hi, H2 [with U{\)' charges Qi and Q2 respec- 
tively; gauge invariance requires Qi + Q2 + Qs = 0] and the /i parameter is 
dynamically generated by the VEV of S as Hs — h{S) — hs/V2. 
The masses of the neutral gauge bosons Z and Z' are 



M; 



Z-Z' 



WMQivl - Q2VI) 



\g[G{Qivl - Q2VI) g[\Qlvl + Qlvl + Qls' 



(2) 



where g, g', g[ are the gauge couphngs of SU{2)l, U{1)y and f/(l)' respec- 



g^ + g'\ 



Vi^2 are the VEVs of with vf 



tively, and 

(246 GeVy. For numerical work we use g'^ ' = (5/3)^' . The Z ^ Z' mixmg 
angle is constrained experimentally to be less than a few times 10^^ (al- 
though larger values are allowed in some cases, e.g., if the Z' has leptophobic 
couplings). We will assume in this paper Qi = Q2 so that the requirement 
of a small Z — Z' mixing will force tan/3 to be close to 1. This constraint is 
less stringent for larger values of Mz'- 

The neutrahno sector has an extra ^/(l)' bino and an additional higgsino 
S (the singlino) besides the four MSSM neutrahnos. The 6x6 mass matrix 
reads (in the basis {B', B, W3, Hi S}): 
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(3) 



where M(, Mi and M2 are the gaugino masses associated with U{1)' , U{1)y 
and SU (2) l respectively. If we assume unification of the gaugino masses at 
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the gauge unification scale, M[, Mi, M2 are in the proportion g'lk[ : ^g'^ : g'^ 
where k[ is a normahzation constant. With this assumption, the neutrahno 
mass matrix depends (for fixed charges) on two unknown mass parameters, 
M[ and Mz' (through fis) and the dimensionless tan/3. In fig. 1 we give a 
contour plot of the lightest neutralino mass in the plane M[ — Mz' for 
tan/5 = 1 and Qi = Q2 = 1- This mass is zero along a line in the half-plane 
M[ > 0. For the study of the region around this line (where x will be quite 
light) it is useful to define the normalized neutralino state A^: 



^ = N 



h 



-=-r^ sin 2(3B' + {Qht' - Qs)S 



+ n{Qs-Q_r^)Hl + r2{Qs + Q.r^)Hl 



(4) 



Z' mixing angle is small due to 



Qic^p + ri,2 = ^^1,2/5 and 



where Q_ = Qic^p — <52'S| (when the Z 
cancellations, Q_ ~ 0); 
This state gives a precise analytic description of the x composition near the 
= line. This is also shown in fig. 1 where the thick solid line delimits 
the area in which the lightest neutralino has a purity P = x ' ^ ^ 0.99. 
In this region Mz' is typically large (M^/ > 300 GeV) so that r is expected 
to be small. To first order in r we can approximate: 



N 



Vr+72 



ricospH^ + sin f3H^) - S 



If we write 



X 



(5) 



(6) 



we find numerically Niiq ~ 0.1 — 0.2, N30 ^ A^4o ~ 0.2 — 0.3 and lA^so] ~ 
0.9 — 0.95 in the region of interest. This is our dark matter candidate and 
it is basically singlino dominated, with a small higgsino doublet component 
and an even smaller B' part. This result is independent of the assumption 
of unification of the gaugino masses and holds as long as they are large 
enough. More precisely, we would have obtained the same eq. (^ provided 
Ml, M2 3> Mz at the electroweak scale. 

The rest of neutralinos (and charginos) have masses controlled by the 
large gaugino masses and fig and are thus much heavier than x- In these 
models the typical soft mass scale is of order Mz' so that the spectrum 
of superpartners is expected to have masses of similar magnitude (a possible 
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Figure 1: Contour plot of the lightest neutralino mass (m^/Gey) in the plane 
M{ — Mz' with Qi = Q2 = 1 and tan/3 = 1. In the region delimited by the thick 
solid lines, the light neutralino has the composition (Q) with a purity > 0.99. 

exception could be the lightest stop because of a sizeable tj?— mixing). The 
fact that electroweak symmetry occurs at a lower scale (and thus Mz -C Mz') 
is a result of accidental cancellations among soft masses (see [Q). The lightest 
scalar Higgs boson, remains also at the electroweak scale (roughly given by 
Mz) while the rest of Higgs states (two more scalars ifg , -f^s , one pseudoscalar 
and a charged pair H^) have heavy masses comparable to Mz'- 
When X is so light that the decay Z XX is kinematically allowed it 
gives an extra contribution to the invisible Z width. The LEP constraint 
5Vinv < 4 MeV is however easily satisfied; the coupling of x to the Z 
boson is proportional to cos 2/5 and this is small, first, because tan/? is 
close to 1 as required to suppress the Z — Z' mixing, and second, because r 
is also small. 

In the region of parameters described above, x is the LSP and, with the 
assumption of conserved i?-parity, becomes a possible candidate for cold dark 
matter. 
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3 Relic Abundance 



The present relic abundance of x's {Q^h"^) is determined by their annihilation 
cross section at freeze-out, which, in the non-relativistic expansion applicable 
in this case, we write as 

^ 2 /^T^ 

(7annV ^ a + -V , (7) 
O 

where v is the relative velocity of x's in the cm. frame. The freeze-out 
temperature Tp (and thus v) can be iteratively computed from the condition 
(obtained from the equality of annihilation and expansion rates): 

_ 0.1Mpi{aannv)m^ 
— = xp - i.n , [fS) 

Ip y/g*Xp 

where Mpi = 1.22 x 10^^ GeV is the Planck mass, g^, the number of relativistic 
degrees of freedom at Tp {^/gl ~ 8 — 9) and 

/ 3 \ 1 

(cTannV) = a + [h - -a] — , (9) 

V 2 / xp 

is the thermally averaged cross section |p. Typically ~ 20 and f ~ 1/3. 
For the range of masses we consider, 5 — 10 < m-^/GeV < 70, Tp is above 
the QCD quark-hadron phase transition and below the electroweak phase 
transition. The neutralino relic density is then 

2_ Px _ 8-77 X lO-i^XFGeV^^ 



X' 



(10) 



where pc is the critical density and h is the Hubble constant in units of 
100 km Mpc~^ sec~^. 

X's in that mass range can annihilate only into standard model fermion- 
antifermion pairsQ (top quarks excluded). The relevant processes are medi- 
ated by Z, Z' and neutral Higgs bosons in the s-channel or by sfermions in 
the t-channel. The cross section formulae for the MSSM can be easily 
generalized to our model. For the particle spectrum and neutralino compo- 
sition described previously, the cross section is dominated by Z' exchange 
with {av = + {bf/6)v^]): 

i2 



2cf 



«/ = "^Pfg't^lQfQxr , (11) 



vr 



^Annihilation into gluons or photon pairs is a one-loop process and can be ignored. 
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and 




(12) 



with 




( 



) 



2 



Am? - Ml, 



\Q)f{2 + if)+2{Qjf{l-ij) 



(13) 



In this formulas, C/ = 1(3) for leptons (quarks) in the final state, ^/ = 



m)/ml, f3f = ^l-if and = QiNl^ + Q2NI + QsNl^ for a x with 



composition as given by eq. (P). The axial and vector U{iy charges of the 
final state fermions are 



Qf = \ [Q'ih) - Q'ifR)] , = ^ [Q'ih) + Q'Ur)] . (14) 



The dependence of cra„„f on these charges introduces some model depen- 
dence in the results. However, annihilation into bb usually dominates so that 
^x^^ depends basically on Q^, which is equal to —Qi if the bottom mass is 
generated by (-f^i). 

Other subdominant annihilation channels are discussed below: 

• The amplitude for Z°-mediated annihilation, which is usually dominant 
for higgsino LSP's, is proportional to {N^q — Nj^), and therefore is doubly 
suppressed in our case: tan /5 ~ 1 implies A^so — A^4o and, furthermore, N^q 
and A^Iq are very small. 

• The t-channel contribution from sfermions, which are expected to have 
masses comparable to Mz', is suppressed by extra powers of small Yukawa 
couplings and/or the smallness of A^i'o, N30, N4Q. 

• The amplitude for XX ^ ^ f f has a gauge part suppressed by 
the smallness of Ni/qN^q, N11QN4Q. There is also a new part proportional to 
hN^Q^N^Q cos P + N4Q sin /5) due to the new Yukawa coupling in the superpo- 
tential (0) and, the sizeable value expected for /i (/i ~ 0.7 from renormal- 
ization group analyses @]), can in principle compensate for the smallness of 
^305 ^40- However, assuming ~ Mz', the contribution is suppressed 
with respect to the Z' contribution by an extra factor m^/vw 

• The same suppression factor appears for xx ~^ H2 — > //. In addition, 
the XXH2 coupling goes like cos 2/3 in the limit uia — rrijjo ^ Mz and thus 
is unimportant for tan/5 ~ 1. 
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• In the limit ^ M^, the third scalar, H^, is singlet dominated and 
does not mediate the annihilation into fermions. 

• More important can be the mediated annihilation; the suppression 
factor m^/vw can be compensated by the fact that m/jO cannot be as heavy 
as Mz' and, in addition, there can be a resonant enhancement of the cross 
section for 2m-^ ~ rrih. Nevertheless, we have checked numerically that this 
channel can also be neglected. First, the CP odd nature of the XX s-wave 
initial state forces to contribute to dannV only to order v^. Second, in 
these models the mass of /i° receives extra contributions and is given by 



m 



ho 



^G^ cos^ 2/5 + sin^ 2/3 + g[^Ql 



V^r + Aradm^O, (15) 



where the last piece, coming from loop corrections can also be sizeable. One 
sees that, even in the case tan/5 = 1, the tree-level mass can easily be as 
large as 170 GeV, and 2my. is never close to the pole for light m^. 

• Finally, as we showed above, the rest of neutralinos and the charginos 
are much heavier than x and "co-annihilation" effects (such as xx' ~^ ff or 
XX^ ff) play no role. 

Figure 2 shows the relic density of x's vs. its mass for different values 
of Mz' from 300 GeV (lower curve) to 800 GeV (upper). For concreteness, 
we have fixed N30 = N^o = 0.25 and N^q = —0.9 in eq. (H) and included 
Higgs subdominant annihilation channels to compute the abundance. We 
use Qi = Q2 = ^ and fix nih = 170 GeV, and = Mz'- In the range of 
rriy, shown, Z' gives the dominant annihilation channel, as explained above, 
so that larger Mz' reduces the annihilation rate making the relic abundance 
grow. For > 70 GeV, /i^-exchange starts to be important and reduces 
significantly the relic density. Smaller values of the f/(l)' charges would also 
tend to increase the abundance but, on the other hand, rrih is lighter in that 
case and will become important for lighter values of m^. 

In the figure, we show only the region Q^^h"^ < 1, conservative limit which 
follows from a lower bound on the age of the Universe of ~ 10 Gyr. Taking 
the uncertainty in the Hubble constant to be 0.4 < h < 1, the cosmologically 
interesting range is 0.01 < ^x^^ ~ ^-^ X'^ form a significant fraction of the 
total dark matter in the Universe. The upper bound f^^^^ ^ ^ ^^^^ ^ lower 
limit on the value of for which x is a good CDM candidate, given fixed 
values of the charges and Mz'- This lower bound increases with Mz' and 
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Figure 2: Relic abundance of x as a function of its mass (in GeV) for Qi = Q2 = 1 
and different values of Mz' from 300 GeV (lower curve) to 800 GeV (upper). 

can be very small for low Mz' without conflicting with experimental bounds 
from decays. 

A non- negligible XX ^ ^' ^ f f ^'^ite usually implies a sizeable Z' pro- 
duction cross section at hadron colliders and one should check that the inter- 
esting region for fl^^"^ conflict with the non-observation of Z' events 
(basically pp ^ Z' ^ e^e~) at the Tevatron CDF and DO experiments. 
The lower limits on Mz' obtained by these experiments are model-dependent 
(all the U{1)' charges and masses of the possible decay products of Z' enter 
the computation) and we do not attempt such detailed analysis. We remark, 
however, that the supersymmetric decays of the Z' (in particular the invisible 
decay Z' — » xX can be very important) reduce significantly the Z' — > e+e~ 
branching ratio, so that the usual limits on Mz' are relaxed and can be 
easily evaded in our model. 

It is illustrative to compare our dark matter candidate with similar light 
higgsino candidates proposed in other supersymmetric models. In the MSSM 
one such light neutralino, with composition x — sin (3H^ + cos (3H2-, was 
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studied in |^ motivated by the SUSY interpretation of the CDF 6677 + fx 
event. The XX^^ couphng, proportional to cos 2/?, is reduced for tan/5 ~ 1. 
This suppresses the contribution to the invisible width and the XX ~^ 
Z —>■ ff cross-section, producing relic abundances in the interesting range. 

In the NMSSM, the minimal model extended by an extra singlet but no 
additional U{1), another light higgsino dark matter candidate was studied in 
[ p!0|] (for more general analyses of NMSSM neutralino dark matter see ||1 1|| ) . 
In this case x — {cos(3H^ + sin/Jif^ + C(S)/N with all three components of 
similar magnitude. The contribution to T^^^, also proportional to cos^2/3 in 
this case, is further reduced by the non-negligible S component. If x — S, the 
relic abundance would be too large because 5"s do not annihilate efficiently. 

Our model is similar to the previous case but with S dominant in the x 
composition. Now, however, x's annihilate through xX ^ Z' ^ ff. The 
rate is not too large due to the smallness of g'lQs and the large mass of the 
Z' boson, resulting in a relic abundance of the right order of magnitude. 

Finally, some comments from the model building point of view are in 
order. The region of parameter space we have examined corresponds to 
M[ ^ Mz' and Mi,M2 ^ Mz- Such hierarchy of masses cannot be easily 
accommodated in models with universal boundary conditions at a high en- 
ergy scale (a GUT scale or the string scale) and would rather point to models 
in which soft breaking is dominated by gaugino masses. Models of this kind 



(in the context of the MSSM) have been considered in the past |jT2[. A renor- 
malization group analysis of the evolution of parameters from low-energy to 
the string scale, with particular attention to symmetry breaking constraints, 
would be required. 



4 Detection Rates 

If x's form the bulk of the dark matter in our galactic halo (with density 
Po ~ 0.3 GeV/cm^ measured by its gravitational effects), we would like 
to estimate the prospects for its detection. We focus on direct detection 
experiments, which are the best suited for small mass WIMPs. In these 
experiments one hopes to detect calorimetrically nuclear recoils in specialized 
materials after elastic scattering with the flux of x's. In the limit N^q = 1, 
the (spin-dependent) x-nucleus scattering proceeds by Z' exchange and the 
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rate of events per day and kg of material is [jTl 

"7.3 X 10^5 



R 



kg ■ day 



(16) 



where is the nucleus mass, {^v)sd takes into account the nuclear form 
factor suppression and velocity distribution of x's, and pl 



C^sd 



16 2Ji-J + '^) 



-m: 



71 



J2 



Y: rf,(Ag(^)(5,) + Ag(")(5„)) 



q=u,t 



(17) 



IS 



In this formula, J is the spin of the nucleus, rrir = 771^171^/(171^ + tti^ 
the reduced mass of the x-iiucleon system, Ag^^^ (Ag*^")) is the quark spin 
content of the proton (neutron), (Sp^n) is the expectation value of the spin 
content of the proton (or neutron) group in the nucleus and finally 



dq = g'l 



2 QqQx 

Ml - 



(18) 



gives the effective neutralino-quark axial coupling. 

For non-zero N^q, N4Q there is also a scalar interaction mediated by 
exchange. The corresponding rate is of the form (0) with a different form 
factor {C,v)sc and asd replaced by 



asc = -ml [Zfp + {A- Z)fn 



(19) 



Z and A are the nuclear charge and atomic number respectively, fp and /„, 
the effective scalar couplings of x to protons and neutrons, are proportional 
to 



Jh 



(P,n) 



iVso (iVso sin (3 + N^^ cos (3) -\ 



m 



2 ' 



(20) 



and TUp^n are the proton and neutron masses, so that fp ^ /„. For heavy 
nuclei, the A"^ factor can compensate for the smallness of N30, N40, and asc 
eventually dominates over the spin-dependent rate. 

We have computed the total rate (spin-dependent plus scalar) for two typ- 
ical materials used in dark matter detectors, '^'^Ge and ^^Si, both with J 7^ 0. 
The results, as a function of tu^, are presented in fig. 3 for the same choice of 
parameters used in fig. 2. In each set of curves, the upper one corresponds to 
Mz' = 300 GeV and the lowest to Mz' = 800 GeV. For Mz' > 400 GeV the 
curves are nearly indistiguishable: the spin-dependent rate is negligible and 
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Figure 3: Detection rates for Ge and Si detectors. 

the scalar rate (independent of Mz') gives all the effect. The small difference 
between the two limiting curves in each set shows that the scalar interac- 
tion provides the dominant effect. The rates obtained are below the current 
experimental sensitivity [|l|. For instance the CDMS experiment, using ger- 
manium, is planning on obtaining a sensitivity of ~ 10~^ event/kg/day after 
a year of exposure. However, the next generation of experiments is expected 
to achieve sensitivities ~ 0.01 event/kg/day and therefore, we conclude that 
x's could well be on the reach of near future experiments. 

5 Conclusions 

In supersymmetric models with an extra U{1) (broken at the TeV scale by 
the VEV of an standard model singlet S) the LSP neutralino can be singlino 
dominated. We have showed that the thermal relic density of this fermion, 
stable if -R-parity is conserved, can be of the right order of magnitude to be 
a good cold dark matter candidate. We have also estimated direct detection 
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rates at typical cryogenic devices in search of halo dark matter and found 
that they are typically small but may be reachable by the next generation of 
experiments. 
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